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The quark-meson coupling model which we have developed previously is ex- 
tended to incorporate the 5 meson. It is then used to study the Nolen-Schiffer 
anomaly and isospin symmetry breaking in nuclear matter. We find that, in com- 
bination with the u-d mass difference, the difference between quark scalar densities 
in protons and neutrons generates an effective neutron-proton mass difference of 
the right magnitude. Finally we find that isospin symmetry breaking in the quark 
condensates can be directly related to the 5 meson field. 
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The Nolen-Schiffer (NS) anomaly [|1], |2[ (sometimes called the Okamoto-Nolen-Schiffer 
anomaly) is a long-standing problem in nuclear physics. The anomaly is the discrepancy 
between experiment and theory for the binding energy differences of mirror nuclei. It is 
conveniently expressed in terms of the quantity AE C : 

AE C = M z> - M z< + 5 np + E c , (1) 

where M z> (z < ) is the atomic mass of the larger (smaller) charge nucleus, 5 np = 0.782 MeV 
is the neutron-proton atomic mass difference and E c is the Coulomb correction associated 
with the additional proton in Z>. Conventional nuclear contributions to the anomaly are 
thought to be at the few per cent level and cannot explain the experimental findings [1- 
4]. The effects of charge symmetry breaking in the nuclear force |5[, especially p-u 
mixing, seem to reduce the discrepancy [4-6]. However, recent investigations of the off-shell 
variation of the p-u mixing amplitude have put this explanation into question [0, § . In an 
alternate approach, using the Nambu-Jona-Lasinio (NJL) model, Henley and Krein |J 
have indicated that the anomaly might be related to the partial restoration of chiral 
symmetry in nuclear matter. More recent theoretical investigations have involved QCD 



sum-rules 10 



In this letter we shall extend the quark-meson coupling (QMC) model |TT| to incor- 
porate an isovector, scalar meson, the 5. We then use it to investigate the NS anomaly 
and isospin symmetry breaking in matter. In this model, which is a natural development 
of the earlier work of Guichon [12 and Yazaki et al. nuclear matter consists of non- 



overlapping nucleon bags bound by the self-consistent exchange of a, u, p and 5 mesons 
in the mean-field approximation (MFA). It provides an excellent description of the prop- 
erties of both nuclear matter and nuclear structure functions |14| . Furthermore, the 



relationship between the QMC model and QHD |L5| has been clarified [[LI]]. 

Let the mean-field values for the a, u (the time component), p (the time component in 
the third direction of isospin) and S (in the third direction of isospin) fields, in uniformly 
distributed nuclear matter with N ^ Z, be a, u, b and 5, respectively. The nucleon is 
described by the static spherical MIT bag in which quarks interact (self-consistently) with 



those mean fields. The Dirac equation for a quark field, if), in a bag is then given by 

[i 7 . d - [rm -V„- 1 -t z V,) - 7 °(K, + ^V,)]^- = 0, (2) 

where V a = g q a, V w = g^u, V p = g q J) and V$ = gjS with the quark-meson coupling 
constants, g%, g q , g q p and g$. The subscripts, i and j, denote the i-th quark in the proton 
or neutron (j=p or n). Here we deal with u and d quarks (i=u or d) only. The bare 
quark mass is denoted by m ; and r 2 is the third Pauli matrix. The normalized, ground 
state for a quark in the nucleon is then given by 

iP ifj (r, t) = Mi/j eM-^t/R,} ^ mXi r/R ) ) ^ (3) 



where 



6i/j = Qi/j + Rj(V w ± T^Vp), for quark (4) 
^7? = 2it*fj 2 (x i )[O i ;(<>. , - 1) + R,m\l2\lx\, (5) 



A/j = v/^'-^O/^ + ^'O, (6) 



with = ^Jxf + (Rjin*) 2 and ^ the quark spinor. The effective quark mass, m*, is 
defined by 

m* = rrii- (V a ± for a ^ quark. (7) 

The linear boundary condition at the bag surface provides the equation for the eigenvalue 

Using the SU(6) spin-flavor nucleon wave function, the nucleon energy is given by 
Ebag + ± \V P for PM, where the bag energy is 

< = E ' n ^ /J ""° + (8) 

with B the bag constant and z a phenomenological parameter (initially introduced to 
account for zero-point motion). Here rii/j is the number of quarks of type i in nucleon 



j. To correct for spurious cm. motion in the bag |16] the mass of the nucleon at rest is 
taken to be 



^ = J(<) a -E"*/jW- (9) 



The effective nucleon mass, M*, in nuclear matter is given by minimizing eq.(|9]) with 
respect to Rj. 

To see the sensitivity of our results to the bag radius of the free nucleon, we choose 
m u = 5 MeV and vary the parameters, B, z and m^, to fit the physical proton and 
neutron masses for several values of the average, free bag radius Ro(= 0.6, 0.8, l.Ofm). 
Since the electromagnetic (EM) self-energies for p and n contribute to the masses we 
adjust the parameters to fit the bare proton mass, M p = 938.272 — 0.63 MeV, and the 
bare neutron mass, M n = 939.566 + 0.13 MeV, where +0.63 MeV and —0.13 MeV are the 
EM self-energies for p and n, respectively fll? |. We then find that B 1 ^ 4 = 187.4, 157.0, 
135.9 MeV, z = 2.043, 1.647, 1.178, and m d = 9.752, 10.01, 10.07 MeV for R = 0.6, 
0.8, 1.0 fm, respectively. The bag radius of the free proton is slightly smaller than that 
of the free neutron, but the difference is negligible. 

For infinite nuclear matter we take the Fermi momenta for protons and neutrons to 
be k Fp and kp n - These are defined by p p = kp p /(3n 3 ) and p n = kp n /(3ii 3 ), where p p and 
p n are the densities of p and n, respectively, and the total baryon density, pb, is given by 
Pp+Pn- The uj field is determined by baryon number conservation, and the p mean-field by 
the difference in proton and neutron densities (p^ below). On the other hand, the scalar 



mean-fields, a and 5, are given by a self-consistency condition (SCC) ||TT| , [I~5|| . Since the 
p field value is given by b = g p pzj \2m 2 p ), where g p = g q and pz = p p — p n , the total energy 
per nucleon, E to t, can be written 

Etot = E [ kFj dkjMf + P + ^ + + ^ PB + (10) 

PM2+) 3 V 2p B 2p B 2ml 8m 2 p p B 

where = 3g^. Then, the SCC for the (= a or 5) field is 
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' 'Mf + P \ 9 <P 



]=p,n 



11 



" (27r) 3 m| 
Using eqs.(|) and (§), we find 

f-^j =-gaX(2C u/p + C d/p )/3 = - g<7 xC;, (12 
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da 



-9a x (C u/n + 2C d/n )/3 = -g a x (13) 



Rn 



{W) R = I x { ~ 2Cu/p + Cd/p) = I x ^ (i4) 

f ' M/; ' ^ ? x (-C u/n + 2C d/n ) = f x Cj, (15) 



where g a = 3g%, g$ = g\ and Ci/j is the scalar density of the i-th quark in j: 




with 



S l/3 - J dr^ lfj - Qi/j{Qi/3 _ 1) + Rjm , /2 - (17) 

Note that while C p has a negative value the others are positive, and that Cd/j is slightly 
greater than C u u. The reason for this is purely quantum mechanical: the u quark is 
lighter, i.e., more relativistic, than the d quark, and hence the small Dirac component 
in the u quark wave function is greater than that in the d quark. In fig.l, we show the 
scalar density factors, Cf, as functions of ps- The dependence of the scalar densities on 
the bag radius is not strong. Using eqs. fll2|) ~ (|15D , the SCC for the <fi field becomes 



(27 



m; 



for ( j meson. (!<■ 



iMf + k 2 _ 

We determine the coupling constants, g 2 . and so as to fit the binding energy (—16 
MeV) and the saturation density (po — 0.17/m~ 3 ) for equilibrium nuclear matter. Fur- 
thermore, we choose gf/47r = 2.82 |l8j and fit the p meson coupling constant so as to 
reproduce the bulk symmetry energy of nuclear matter, 33.2 MeV. We then find that 
gl/Air = 19.4, 20.8, 21.1, g 2 jAx = 1.69, 1.26, 1.08, and g 2 /4rr = 5.51, 5.74, 5.78 for R 
= 0.6, 0.8, 1.0 fm, respectively. (Here we take m a = 550 MeV, m w = 783 MeV, m p = 
770 MeV and mj = 983 MeV.) The present model gives a good value for the nuclear 
compressibility - around 220 MeV. 



Since we have seen the a meson field in Ref. [11], here we show only the mean field 



values of the 5 meson in fig. 2. We define the proton fraction, f p , as p v j pn- Clearly the 6 



meson field has a symmetry property, 5(f p ) = —8(1 — f p ) in the case where m u = rrid, so 
that in this case 5 = for f p = 0.5. However, once the quark mass difference is introduced 
this symmetry is lost; even for f p = 0.5 the 5 field does not vanish but has a negative 
value. It is easy to see why it is negative: from eq . flUf) , the a and 5 fields at low ps can 
be expressed as 

g a a~ (^f (C; Pp + C° nPn ), (19) 
g s S^~(j£j\c s p p p + C s nPn ). (20) 

Therefore, when f p = 0.5, g$5 = —\ (^) (C p + C^)pb which is less than zero because, 
as explained below eq . (|17D , C p + > 0. 

At this stage we are ready to consider the NS anomaly. As the first-order approxi- 
mation, the binding energy difference, eq. ([]]), would be treated in terms of the difference 
between the effective masses of p and n in matter: AE C = A° np - (M* - M*), where A° np 
= 1.293 MeV is the n-p mass difference in free space. Since the effective masses, M* and 
M*, are functions of ps in our model, we may use the local density approximation (LDA) 
to calculate AE C for finite nuclei. For medium mass nuclei, AE C has been measured us- 
ing mirror nuclei ||. On the other hand, for large nuclei (i.e., neutron rich nuclei), the 
experiments have been performed by measuring analog resonances in (p, p') and (p, n) 
reactions ||19| . In LDA the effective n-p mass difference in the nucleus A, A^ p , A , may be 
calculated as 



A* 



JdrA* np (p,f p )P(r), (21) 



where p describes the core density of the mirror nuclei or the nuclear density for large 
mass nuclei, which is taken to have a Woods-Saxon form: po/(l + expf 11 ^]) with p the 
normalization factor and a ~ 0.54/m. Here we take R so as to reproduce the rms radius 
of the core nucleus for medium mass nuclei || or R ~ l.l2(fm)A 1 ^ - 0.86(fm)A- 1 / 3 
for neutron rich nuclei. Furthermore, we choose A* np = M*(p, f p ) — M*(p, f p ) with f p = 
0.5(Z/A) for mirror (large) nuclei. The probability distribution, P(r), of the last nucleon 
in mirror nuclei is given by the harmonic oscillator model with the radius parameter 
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chosen to fit the rms radius of the valence nucleon p0 |. (For large nuclei we use the 



nuclear density distribution itself for P(r).) We show the bag radius dependence of A* p 
in fig.3. It gradually decreases as ps grows, becoming negative for high densities, and 
it does not much depend on R Q . We can easily see why it decreases: since the effective 
nucleon mass can be expanded at low densities as 



Mj^Mj-grCp + ^Cfa (22) 



we find that A* p = A° np - g a {C° n -C°)a+ \g 5 {C 5 n - C 5 p )5 when f p = 0.5. Because Q > C£ 
and 5 is numerically much smaller than a, the n-p mass difference decreases as ps goes 
up. The numerical results for AE C (= A° p — A^ p , A ) are displayed in table 1. The present 
model provides a good overall description of the variation of AE C with mass number. 
We can see the noticeable shell effects arising from the differences in the average density 
felt by the valence nucleon. The calculated values are of the right order of magnitude to 
resolve the anomaly. Note, however, that a complete treatment of the anomaly, including 
the spin-orbit effects, goes beyond the local density approach used here. 

Next we wish to use the model to study the quark condensates in matter. Within 
the QCD sum-rule approach these play a very important role in a wide range of nuclear 
phenomena. The difference between the z-th quark condensates in matter, Qi(pb), and 
in vacuum, Qi(0), can be related to the meson fields through the Hellmann-Feynman 
theorem [[11], 0: 

d£ dM* 

Qiips) = Qi{pB) - Qi{0) = J2 ^ttt _ 3 + (meson contributions), 

j^n dM j dm i 

d£ d£ 

= n i/ P Ci/p-^; + n i/n C i/n -^^ + 0(p 2 B ), (23) 

where £ = psEtot and the meson fields contribute 0(p 2 B ). We take up the terms of O(ps) 
to see the quark condensate at low density. Defining the average of the quark condensates, 
2) by \{Q U + Qd) and their difference, SQ, by Q<i — Q u , they can be expressed as 

Q(Pb) * \ (^) {g.o) * \{C; Pp + G>„), (24) 



and 

/ \ 2 

1115 \ („.x\ ~ n8„ 1 n& 



SQ(pb) — -2 [~J (g S S)^C d pPp + C A nPn , (25) 



where we used eqs.(^) ~ (pO|). It is worth noting that at low density Q and 5Q can 
be described by the a and 5 fields (respectively) alone. Since we have already discussed 
Q in Ref.|llJ, we shall concentrate on 5Q. In particular it can be related to 7, which is 



a measure of isospin symmetry breaking in the quark condensate in matter, defined by 
7 = Q(i(pb)/Qu(pb) — 1- Using eqs.(|20|) and (|25D , we find that at low density 7 is 



2 

7 " 70 + & " 70 + mm (if) " 70 - ift^ (c ' p ' + c " p " ) ' (26) 



where |Q«(0)| ~ (225MeV) 3 and 70 is the value of 7 in vacuum: 70 = —(6 ~ 
9) x 1(T 3 |22[. Therefore, one finds that 7 ~ 7o - 1.0(1.3) x 10- 3 (p B /p ) for R = 
0.6(1.0)/m when f p = 0.5. Since 7 is directly related to 5, 7 depends not only on ps 
but also on f p , which is quite different from the result in the NJL model ||10|| . Our result, 
however, seems natural because the origin of isospin symmetry breaking is very clear in 
our model: the 5 meson couples differently to the u and d quarks and the difference of 
their effective masses in matter is completely controlled by 5 (see eq.([7|)). Note that 7 
can be rewritten in terms of the quark effective mass: 7 ~ 70 + |q 2 (o)| (^f ) 2 (^ m d ~ ^ m u)^ 
where 5m* = m\ — mi. 

Finally, we shall consider a potential experimental signature of isospin symmetry 
breaking. Since the u-d quark mass difference leads a difference between small Dirac 
components in u and d quark wave functions, the axial-vector neutral-current coupling 
constants of p and n, g p A and g% seem interesting f23| . Using eq.(|2|), the ratio of g\ to g p A 



is given by g n A j g v A = (4G d / n + G u/n ) / ' (AG u/p + G d/p ), where G i/:j = M?^ f drr 2 {j$(xir/ Rj) - 
iPi/jfifar/Rj)] = {2Ql J +AQ l/j R J m* j -3R j m*)/3{2Q i/j {Q i/j -l) + R :j m*). In symmetric 
nuclear matter, we find that g A /g A = 1 + fr/ with 1] ~ 2(4) x 10 _3 +0.2(6) x 10~ 4 (p B /p ) for 
Rq = 0.6(1.0) fm. The dependence of rj on p B is quite weak. At the present experimental 
level it is not possible to see it, but one may hope in the future. 

In summary, we have applied the QMC model involving a, u, p and S mesons (all 
mesons whose masses lie below 1 GeV) to investigate the NS anomaly for finite nuclei. 

8 



We have also used it to investigate isospin symmetry breaking in the quark condensate in 
nuclear matter. The result is a very reasonable description of both the anomaly and 7 in 
nuclear matter. Their physical origin is very clear and the differences between the quark 
scalar densities in the proton and neutron, which are generated by the u-d quark mass 
difference, are essential to understand them. 

This work was supported by the Australian Research Council. 
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Figure captions 

Fig.l Scalar density factors, Cf, in symmetric nuclear matter (Rq = 0.8/ra). The solid, 
dotted, dashed and dot-dashed curves are for C£, C°, \C^\ and C*, respectively. 
The values of the coupling constants are given in the text. 

Fig. 2 Mean field values of the S meson for Rq = 0.8/m. The dotted curves are for 
the case where m u = m d = 10 MeV. The solid curves display the full calculation 
involving the quark mass difference. The top (middle) [bottom] two curves are for 
f p = 0.7 (0.5) [0.3]. 

Fig. 3 The n-p mass difference in matter (f p = 0.5). The dotted, solid and dashed curves 
are for the calculations with R = 0.6, 0.8 and 1.0 fm, respectively. 
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Table 1: AE C for several finite nuclei. Energies are quoted in 



R (fm) 


0.6 


0.8 


1.0 


observed discrepancy 


15 _17 N 


0.29 


0.32 


0.34 


0.16 ± 0.04 


17 F _17 


0.22 


0.25 


0.27 


0.31 ± 0.06 


39 Ca- 39 K 


0.36 


0.41 


0.44 


0.22 ± 0.08 


41 Sc- 41 Ca 


0.34 


0.38 


0.41 


0.59 ± 0.10 


120 Sn 


0.72 


0.83 


0.87 




208 Pb 


0.78 


0.91 


0.95 


~ 0.9 
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